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ABSTRACT 



A probabilistic model is developed for treating the 
problem of optimizing the intended trajectory and associated 
height of burst of a missile-warhead. The angle of re-entry 
is treated as a control variable. Computational techniques 
which may be operationally acceptable are described. 
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I INTRODUCTION 



Consider the problem of setting the fusing mechanism of a 
warhead aboard a ballistic missile in order to optimize the 
probability of destroying a target. 

The general problem as stated above becomes extremely 
complicated when one considers such things as missile in-flight 
failures, warhead failures, back-up fusing mechanisms employed, 
type of warhead and type of burst, distance from launcher to 
target (which will affect range and Azimuth errors), geodesic 
error, (i.e., to what degree of accuracy do we know the 
distance to target.), and more frustrating, the problem of 
partial kills, i.e., the partial destruction of the target or 
a portion of the target for a period of time less than the 
duration of hostilities. 

In order to reduce the problem to something that can be 
handled, it will be convenient to cover briefly three topics: 
the concept of ballistic coefficient, the reflected (mach) 
wave phenomenon and the definite range law (cookie cutter) . 
Ballistic Coefficient : 

For our purposes, we may think of the ballistic coef- 
ficient as a parameter which describes how much a re-entry 
vehicle will be affected by the atmosphere (see Fig. 1). It 
is proportional to the ratio of the weight of the vehicle to 
its area, i.e., the square of the maximum diameter (see Ref. 4); 
e.g., for a given area, an increase in weight will tend to 
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"straighten" the trajectory whereas a decrease in weight will 
result in a slower re-entry and thereby increase the time spent 
in the gravitational field. The result will be a more curved 
trajectory. 

Mach Effect : 

When a warhead is burst in the air close to ground (e.g., 
300 feet for a 1 KT warhead), there is a region of unusually 
high overpressures due to the merging of the incident and re- 
flected waves (see Ref. 1 and Fig 2). This region is known as 
the Mach Region. Because of this the curves in Fig. 2 are not 
monotonic but have the unusual shape exhibited there. 

Definite Range Law : 

This is the "cookie-cutter" concept that is frequently 
used. Because of its simplicity it seems to give useful 
results for many objectives. 

A couple of examples should be sufficient for our purposes. 
Consider a point target in space at which we are firing a pro- 
jectile with a kill range of one mile. Imagine the target to 
be at the center of a sphere of radius one mile. If we detonate 
this weapon anywhere within the sphere we consider the target 
destroyed; anywhere outside the sphere is considered a miss. Of 
course, no weapon behaves in such a manner. "Partial kills" are 
generally the rule rather than the exception, and they are very 
difficult to quantify. The definite range law is a mathematical 
convenience. Presumably, this range can be regarded so that a 
controlled fraction of partial kills will be treated as misses 
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Figure 3.67&. Peak' overpressures on the ground for! a 1-kiloton burst (high*pressure range). 



FIGURE 2 



“FJg ,~2^epYesents a destruction -curve for various overpressures,' .of a ' 'V 
1 KT weapon. For a weapon of different* yiel4 M 'the'/ o’rdi'tiete, and ’abscissa \ 
are each multiplied by W J .. See Ref. 1. For exalnple^’a 60,0, foot high burst 
of a 1 KT weapon, the 20 P.S.I. curve ( will 1 , extend to about 940 ^ fee t' from \ 
ground zero. Similarly, a 1 MT weapon detonated* at» j'6 ,000' 5 feet will have a I? 
20 P.S.I. overpressure to. about 9,400, feet. * '* * ! ' ' 



TECHNICAL ASPECTS OF BLAST WAVE PHENOMENA 



(and a controlled fraction will be scored as hits). 

As another example consider Fig. 5. If we imagine a target 
at (0,0) and estimate that it takes about 20 P.S.I. of over- 
pressure for destruction, the 20 P.S.I. curve shown could be 
considered as a "cooker-cutter"; i.e., a burst just inside 
will destroy the target and one just outside will not. 

In light of the above discussion, the following three 
simplifying assumptions will serve as the basis of this paper: 

1. The re-entry vehicle has a sufficiently high ballistic 
coefficient so that the trajectory may be assumed to be a 
straight line in the neighborhood of the target. Moreover, 
the vehicle is launched from a sufficiently great distance so 
that all possible paths may be assumed to be parallel in the 
neighborhood of the target. 

2. The overpressure contours of Ref. 1 (and hence of Fig. 2) 
are valid and the target damage is a function of overpressure 
alone . 

3. A definite range law holds and 20 P.S.I. is necessary for 
destruction. 

The general objective is to select an optimum height of 
burst. Further examination of the problem reveals that some 
other parameters enter into the picture. The angle that the 
vehicle f s trajectory makes with the surface can be selected 
in advance. There is an advantage to small angles since then 
the missiles are more difficult to detect (and hence defend 
against). Larger angles generally lead to higher kill 
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FIGURE 3 



DESTRUCTION ENVELOPE 



probabilities. The capability to vary this parameter will make 
defense more difficult. Given a trajectory angle, the optimum 
burst point becomes a function of height of burst and aim point 
(intersection of the intended trajectory with the surface). 
Locating the target at (0,0,0) and using the Cartesian system 
wi th 



)C,= downrange distance 
X,= crossrange distance 
X3 = vertical 
and introducing parameters 

angle the trajectory makes with the surface 
? = distance from aim point to target 
= intended height of burst 
(see Fig. 3), the problem now may be stated: 

Given-0" , choose the pair ( ^ ,^/A ) so that the probability 
of target destruction is maximized. 

This paper must be viewed as a pilot study for solving the 
problem outlined above. It will be noted that a large amount of 
detailed work remains to be done in order to provide good answers. 
The methods herein can serve as outlines and guides. 

One of the more important results is the development of a 
graphical computation system using dividers, graphs of over- 
pressure functions and correction curves. Thus, the optimum 
solution may be obtained under operational conditions without 
the need of a digital computer. 
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The organization of the report is as follows: The mathe- 

matical model and nature of the problems of mathematical 
analysis are presented in section II. The computational problems 
are described in section III. The operational hand approxi- 
mations and corrective curves are given in section IV. The 
conclusions appear in section V. The Monte Carlo technique 
employed, along with the computer program, appear in Appendices 
A and B respectively. Appendix C contains a tabulation of 
the parameter values determined by the graphical method de- 
veloped in section IV. 
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II 



CONSTRUCTION OF THE PROBABILISTIC MODEL; 

MATHEMATICAL CHARACTER OF THE PROBLEM 

We have already introduced the coordinate system; the 
angle the trajectory makes with the surface; the intended 
height of bursty^A ; the aim point ( , 0, 0). The trajectory 

is actually a random phenomenon and by assumption 1 can be 
characterized by the impact point and the incidence angle . 
The actual height of burst is also subject to random errors. 
Thus, we let 

(Yi , Y£, 0 ) be the random vector representing the impact 

point , 

(Xi , X 2 , X 3 ) be the random vector representing the burst 
point 

A sample situation is given in Fig. 4. 

It is assumed that: 

1 . Yi,Y 2 has a circular normal distribution with mean vector 
^ y 0 and covariance matrix 0 ^ I where I is the identity 

matrix. 

o 

2. X 3 has a normal distribution with meany.,*, , variance^ and is 
independent of Yi, and 

From Fig. 4 it is seen that Xi = Yl - X 3 cot 

X2 = Y2 
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It follows that the probability density function of the 



detonation point is 






The ballistic error parameter 07 , and the fusing error parameter 



CT'should be determined from independent experiments. There are 
many types of warheads and fuses and we merely carry (T , and (T 7 
as "known" parameters. 

Since 20 P.S.I. will destroy the target, the event of 
target destruction can be visualized in the , X 2 , X 3 ) space 
by rotating the curve in Fig. 2 about the vertical axis. If 
we call this set A, then the probability of target destruction 
may be expressed as the triple integral of the density function 
(2.1) over the set A 



It may be well to clarify a small logical difficulty. The 
probabilistic model permits negative values for the height of 
burst. These are physically impossible and should be interpreted 
as a duded warhead that impacts the surface without exploding. 

The optimum solution can be obtained by taking the appro- 
priate partial derivatives of ( 2 . 2 ), equating them to zero and 
solving the resulting system of equations. Thus 



( 2 . 2 ) 





(2.3) 




(2.4) 



bP 

c )/>v 



= O 



£jL 
t d§ 



- 0 



Examination of the character of equations (2.4) is facili- 
tated if we make the change of variables 



( 2 .5) y, -- 

<T , 



1 



V - '^' L / ~ ^3 

y L- —r ) h ' ~=L 

(T. ' J r 



and let 






£ 

<r 
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In this coordinate system the set A is transformed into a set B 
the exact nature of which is difficult to visualize (see Fig. 2), 
but we can say B is bounded since the transformation is linear. 



In the new system, equations (2.4) may be expressed as 

Cry. -fA Yt f(y, -/*>'] 



\ ^ff j 



- 0 



(2.6) y 

(-lY'- 

u# 
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jff(yrr) e 1 

-frowsy 

•y.-f) e ' 

Geome trically , the problem now may be viewed as follows: Given 

a set B, choose a centering point ( , O ) so that the 

probability of B is maximized. 

All solutions of (2.6) are critical points which may be 

r* 

local maxima, minimum or saddle points. The nature of these 
points may be examined by means of second order derivatives. 

letting ,,x -XiYrdWm-vy 

< 2 . 7 ) MlM (iir) £ 



we can write i 

( 2 . 8 ) ^ -- fff L ( W '- 0 4 (\ 4 , y, i H < * 

(2-9) -- fff t'f'A'i 

(2.10)0,5 "fff[(Y,- f )‘-lJ J 
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A critical point will be a local maximum if 



(2.11) D D 



■ % >° 



and 



(Ref. 3 p. 232) 

It can generally be said that D^y (and D r y , for the same 
reason) is negative since the variance J j f K ( t 



of Y 3 is one and ( 2 . 8 ) compares a truncated (by B) version of 
this variance with unity. The remaining part of (2.11) is a 
difficult question, however, and may require more specific in- 
formation about the set B. 

The set B plays a role in determining the number of critical 
points. For example, if B were a sufficiently long dumbell 
shaped region, there would be at least two solutions to ( 2 . 6 ). 

It would be interesting to examine the question of number of 
critical points if B were convex. With this condition it may be 
possible to show there is only one. This contingency was not 
examined since the set A (and hence B) is not convex. 

Since the set A is formed by rotation of Fig. 2 about the 
X 3 axis, it can be shown that 



( 2 . 12 ) S { -o 

A 

The same must be true in the transformed coordinate system, and 
it follows from ( 2 . 6 ) and ( 2 . 12 ) that (£ , O^U) will be the 
center of gravity of the conditional probability density given 
(Y]^,Y 2 ,Y 3 ) belong to B (if f and ^ are unique solutions of 



( 2 . 6 )). 
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COMPUTATIONAL PROBLEMS 



Computer techniques are required for the solution of the 
system of equations (2,4) , the examination of the critical 
points ( 2 , 11 ), and the determination of the optimum probability 
of target destruction (2.2). It was decided to use Monte Carlo 
techniques (see Appendix A) in this study, and the first step is 
to characterize the set A. The curve in Fig. 2 may be viewed 
as giving distance as a function. 

This curve is not given in analytic form. Also, note of 
height, i.e., g(X^ ) that there is a small interval of X 3 
values for which g(X^ ) is double valued. It was decided to 
ignore the double valued feature and fit a polynomial of 
order 7 by the method of least squares, using 23 values read 
directly from the graph. This fit appears in Fig. 8 . Using 
the polynomial as g(X 3 ) , the set A may be represented by 



The integrals in (2.2), (2.4) and (2.11) all are similar 
in nature. A rapidly converging iterative technique may 
generally be available for the solution of the system (2.3). 



(3.1) 




2 2 
Xi+ x 2 £ 




Letting 




the equations (2.4) may be written 





Let functions f ) and g2(/ A / / ) be defined so that the above 

equations are 

, f-fai/.f) 

Then, using an initial approximation ofyK^ ^ we can use 
the iterative scheme 

(3.3) /s.. -- fa (AS £, ) _ /, + , fa (/S L ) 

provided it converges. 

The examination of the question of convergence usually 
goes along the following lines: If the series 



(3.4) 






'll 






converges absolutely, ^hin and ^ converge to finite 
numbers which will be the solutions of 2.4. This will take 
place if a constant r ( 0<r<l) can be found such that 

*A+i ~AI - i~(l/ A /-/ u ,-,i + Ifj-ft-'l) 

°' 5) I Vi - 6 I - AA'/V'I + I 

Tne question of convergence of the scheme is more easily 
examined after making the transformation (2.5). Let this be 



for then 



and 



done and write ~ 

<? 

tt / 7 T Yj('',-f-)fi('v0(y r /i*) 



and then 
(3.6) 
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Proceeding as indicated, 

o.r> \f» h rh\- I ^ r(J~>30Cn- f 

¥ \(t~ y >i(Y.-r,)0(Yi)4(' / r/ A i ll 

+ l£. f/f r i$(Y.-fj)0(Y>)0(Y, Yrftf ,- I 

"Ifr i.if/f Ys0a-!,)<lM<0( y rA-') 

* fit {/ f I '/, I (Z5 fy , ) 01Y, -A- JIM - fj 

From the theorem of the mean 

0f7v- /fry - /rn ***4'{t) 

it is easily shown that ( 

< 3 - 8 > M'W/i ^C' Ts (°</ 

and this value may be used in 3.7. Thus we can write 

| 4 (fr/'il- 0 (/’-/*,-,) I s P/A/V'l 

1 0(t- f ,)- (/(Y.-f,.,)! iC\^ - f,.,l 

Sharper bounds may be available if the maximum is 3.8 and the 



are constrained to be in the set B. The 



points ( f ,o ,fl) 
second term of 3.7 can be estimated in a similar fashion: 

I ff/m- f r ,) p(Y)/>(Y s '/T / (Y,-r,$(M0(Yr/ L '/)]l 



which can be treated as above. 
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Collecting, we have , r . , 

+ 1 , p i f, - f,-i r/( f/,; ■ /V, J 

(3.10) 

♦ftf/ | %\$w,-f,)0('L)p(y, ■/<,«)] ^ f Pl/'r/'r' 1 
[H? (Y, ( Y*J + PI f, - f ; ., I ffftf (yjffty ■/*, )J 

which will have the form (3.5) if it can be shown that the 
coefficients are less than one. This is a difficult analytic 
problem and is beyond the scope of this paper. The technique 
(3.3) was used in our computations and lend to answers that 
were stable in the light of our capability to compute by Monte 
Carlo. The signs of the determinants (2.11) were computed at 
the solution values and all were found to be positive. 
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FIGURE 4 

SKETCH OF TRAJECTORY AND CO-ORDINATE SYSTEM 
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IV 



A SIMPLIFIED APPROACH TO OBTAIN AN APPROXIMATE SOLUTION; 
GRAPHICAL METHODS AND CORRECTION 



If the down range and cross range errors are assumed to 
be small, the problem can be reduced to a simple two-dimen- 
sional model since we intuitively feel that a 'best 1 trajectory 
is one that passes directly over target. 

At this point we will assume we follow our intended 
trajectory with probability 1 and find values of JA and > that 
optimize this conditional probability of destruction. 

If Fig. 2 were reflected about a vertical axis, Fig. 5 
will result. Again, assume the target to be at (0,0) and we 
wish to detonate our weapon inside the 20 P.S.I. envelope. 

We will make the following assumption to supplement those 
in section III: 

The warhead passes directly over target in a straight 

line trajectory making an incidence angle of0-with the surface. 

The collection of possible trajectories will intersect the 20 P.S.I. 

envelope as shown in Fig. 5. Denote the lower intercept (of a 

given trajectory) by a and the upper by o< . 

The maximization takes place in two stages. First, holding 

■0 and T fixed, find the best value ofyL{ as a function of f . 

Second, vary \ until the maximum is obtained. To implement this 

procedure, let . , 

, r-J 

(4.1) F ( M) ~ "- 1 - p 

r m<r J* ^ °v 
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i tjsjasi kis/kU Abr&uio ttu/vai wa.vju 




SAMPLE OF METHOD EMPLOYED IN ISOLATING OPTIMUM-INTENDED BURST POINT 



or 



Taking the derivative and equating to zero 

<*■’> />>- C«Tr 

Making the substitution t = Y-jj 

, r 

r±£ tl 

<r 

- <v 4 U =* 

r -t (¥; 1 
€ - € 



-(7 



Since a function and its log have a miximum at the same points 






Simplifying and solving for jj 



z _ ^ 



/A - 

'2 ( - <5c) 



(4.3) 



/*« 



4 / 



- + a 



o( 7 ^ 



(4.4) 



The problem is degenerate for the case o( = a since the 
probability of detonation is zero. It is interesting to note 
that the above solution is independent of (f and the velocity 
of re-entry. 

To verify we have a maximum, consider 

, , . iLJf w -(M)' 

1 ~ietr (fflz ^^/ + Lwrr^ Xr <*y 
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Making the substitution t = Y - p and simplifying we obtain 



after one integration by parts we see 

- (JL2> s ( £ 
urn? z 

which is true for all ^ } Q 

The value of F at the point ^imax which is the probability 
of destruction given a particular line is given by 

iY ' *%/ 

Letting Tt--- y- 

f w - r® s €■ 



S l 
^ <r 



or 



(4.5) 1= (*&) S' ^ 



Q (-go 

t <r 



- oo 



UTrr 



» X 

<? ^ At - 1.0 



So far we have the optimum height of the IBP, but still 
must find out where to plan the detonation, i.e., to select f . 
Since we are interested in maximizing the probability of 
destruction (4.5), we see this is clearly equivalent to 
maximizing ( d -a) since the function is monotonic. Maximizing 
( d -a) is in turn equivalent to finding the longest line that 
can be drawn through Fig. 5 at an angle & . 

In Fig. 5, a series of parallel lines were drawn roughly 
every 20 feet at an angle-© to the horizontal and the longest 
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line determined by the use of dividers, and choosing the one 
with max ^ a max i mum occurred between lines it 

was deemed to be at the mid-point of the two which gives an 
"accuracy" of 10 feet for a 1 KT weapon or 100 feet for a 1 MT 
weapon . 

By assuming the measurements were accurate in the sense 
that one can differentiate between a line being longer than an 
adjacent line (which is reasonable since the vast majority 
were easily discernable) we can say with probability one, the 

V} 

maximum error is 10 W feet for a W KT weapon. 

Measurements were taken in this manner for several values 
of and the resulting values of ^Imax were plotted against . 
Readings were concentrated at points where the curve appeared 
to be rapidly changing. A similar procedure was used for values 
of 6* vs / and ^ vs J s . in order to interpolate polynomials 

of degree 1, 2, , 20 were fit in a least- squares sense to 

the data points. The degree chosen from the one with the 
smallest ^ I Y A ~ G A j J where = emperical value 

and Oi = computed value. The above criterion was chosen since 

it was planned to interpolate values of "G~=5, 6, 7, , 90 

and the measurements were felt to be fairly accurate to start 
with and large deviations in the fit would only tend to 
aggrevate the original data points. 

The polynomials were fed into the computer and evaluated 
at the points &= 5, 6, 90. The results were tabulated 
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in table 1 where: 



jimax = planned burst point measured vertically along HOB axis 
XCOORD = planned burst point measured horizontally along 
"DISTANCE FROM GROUND ZERO" AXIS 

= Angle trajectory intersects ground measured from the 
horizontal 

f = Aim point measured positively to the right along "DISTANCE 
FROM GROUND ZERO" AXIS 

REMARKS CONCERNING THE GRAPHICAL APPROXIMATION 

The most interesting discovery is that the 'best* IBP is 
never over target; for low angles of re-entry, the burst is 
planned prior to target; for steep angles of re-entry, the IBP 
is past target. 

The sharp discontinuity in XCOORD at about 28° is caused by 
entry into the envelope in the region (850, 700) of Fig. 5 where 
the slope of the envelope and the slope of the trajectory are 
nearly equal . 

We will next apply the correction discussed in section III. 
The method developed above was used as the first approximation 
for the iterative technique developed in section III. The 
details follow. 
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Let 



- iT_f ^ 



2<r, 






?- f-f 

Pf -' |f Jr x (Arf^X 3 corer ^'0 



( a ,-/-; 



Ut f. * |A,f , F;--Af 

E, - f Pf cot&E s =o 

E r/ uP ;0 



Solving for ytA and y 

- £ 

P 



yU r _cj 



p. £, +/ApC0T& 
' J ~ p 



let yu | ^ f i 



= First approximation obtained by using the 
method of section IV 



- y i fin ‘ 

In general, about 10 iteratives were needed for# > 10° 

and about 20 for# = 10°. 

Data was collected for *# = 10°, 20°, 80° for two 

sets of the standard deviation or; *7 ) . Table 1 contains the 
percent of error in^M and J and Fig. contains a plot of # vs 
percent of error. 

Thus, operational procedures can be developed that do not 
require a high speed computer. The first approximation is 
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0 ~ = 500 



< r \ = 200 C = 200 , = 500 



-e- 


A 


f 


J* 


f 


10 


+ 20.8 


+ 0.1 


+50 


+25 


20 


+ 7.5 


+ 1.34 


+ 48.9 


+71 


30 


+ 4.28 


- 1.8 


+ 12.5 


+ 5.9 


40 


+ 7.9 


- .134 


+ 3.54 


- 24.2 


50 


+ 1.5 


- 4.8 


- 6.15 


-45 


60 


+ .065 


- 2.05 


- 5.6 


- 41.5 


70 


+ .43 


- 1.9 


- 7.1 


-65 


80 


+ 1.3 


- 3.2 


- 5.8 


- 14.75 



TABLE 1 

PERCENT OF ERROR IN yV AND f BY USING FIRST APPROXIMATION 
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<r = 200 
<r, = 500 



-e- 


M 


f 


PKILL 


10 


453 


2504 


.419 


20 


431 


1146 


.582 


30 


432 


719 


.652 


40 


440 


508 


.684 


50 


438 


369 


.692 


60 


438 


256 


.699 


70 


435 


150 


.712 


80 


438 


84 


.716 



< r ~ = 500 

<rr = 200 





/ 


? 


PKILL 


10 


365 


2002 


.227 


20 


315 


769 


.402 


30 


391 


658 


.499 


40 


459 


661 


.569 


50 


472 


538 


.615 


60 


468 


429 


.637 


70 


463 


306 


.642 


80 


471 


155 


.649 



TABLE 2 

PROBABILITY OF DESTRUCTION AS A FUNCTION OF RE-ENTRY 
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GRAPH OF PERCENTAGE ERROR IN PREDICTING PROBABILITY 
OF DESTRUCTION VS ANGLE OF RE-ENTRY 



obtained by the use of dividers applied to graphs like Fig. 5. 
The results are corrected by curves like Fig. 6. 

As was expected the method of section III gives a more 
accurate first approximation when the ballistic sigma is small 
Probabilities of destruction can be predicted more accurately 
for higher values of In general, the method of section IV 

gives values of yWand f* that are low and nearer to target for 
small -g*- ( 0- < 40°) ; high and further from target for higher 

This has high intuitive appeal since if re-entry is at a 
small angle we tend to shoot for the point in Fig. 5 where the 
destructive envelope intersects the ground. But if ballistic 
errors are allowed, we run a high risk of having the vehicle 
impact with the earth before ever reaching target. Similarly, 
for steep angles of re-entry we tend to aim our errorless 
vehicle to pass thru the highest point on the envelope since 
the vertical distance traversed is a maximum. But if we allow 
errors, we run a high risk of passing to the left of the 
envelope of Fig. 5. 

Appendix C contains a list of starting points </•. > 

for^* =5, 6, 7, , 90 along with associated destruction 

probabilities. If one chooses to use the curves given in 
Fig. 5, the destruction probabilities given may serve as a 
check on the final probabilities, the percent of difference 
decreasing as (J*j decreases. 
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Outside of Appendix C, all probabilities were computed 
using the polynomial fit to the curve since the curves were 
not believed to be accurate in the first place and we were 
interested primarily in technique. 

Graphs of the polynomial used and the actual curves are 
found in Fig. 7 . 
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CONCLUSIONS 



The optimum probabilities of target destruction, and the 
optimum values ofyl^and J~appear as a function of & in 
Table 2 for parameter (<T) ) =(200, 500) and (500, 200) ft. 

Notice that the fusing error parameter (P is the more 
important parameter in terms of increasing the probabilities. 

This paper should be viewed as a pilot study. The method 
should be sharpened in the following ways: 

1. The overpressure envelope curves (e.g., Fig 2), should be 
established as a function of yield and mathematical methods of 
characterizing them (e.g., Fig 7), should be improved. 

2. Sharper numerical integration techniques should be applied. 

The present Monte Carlo method has a 3 standard deviation 
tolerance of + .015 for the probability of destruction. 

3. The technique should be applied to a large number of yields 
using realistic delivery errors and fusing parameters ( (Tj and (T ) . 

4. The hand method for operational computations needs proper 
development, particularly in the area of sharpening the 
correction curves (e.g., Fig 6 and Table 1). 

For the parameter values used in this study, the probability 
of destruction was not very sensitive to changes in the values 
of ^ and ^ . This could change with other parameters. For 
example, a 5 MT weapon delivered from 2000 miles 
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may well have a delivery error of say (T| = 5 miles. If 50 P.S.I. 
is required for target destruction, it may be necessary to 
determinedly and ^ very precisely. 

It was felt the assumption of a straight line trajectory 
is fairly realistic and that every moderately large deviation 
from a straight line will yield approximately the same results 
since the major dependency is simply in the fusing CP which is 
largely dependent on the particular warhead and not on the 
trajectory. Since curved trajectories mean lower speeds, the 
fusing error ^"’"will be smaller. It follows that the probability 
of destruction will be higher. Thus, the use of straight line 
trajectories will yield probabilities which may be used as 
lower bounds. 

The normalicy assumption in the fusing (p is felt to be 
pretty good if the warhead is fused by a timing mechanism or 
a radar, but barometric devices are notoriously non-normal. 

The circular normal distribution assumed for the trajectory 
is probably not realistic but is of little consequence since the 
model is broad enough to include different cross range and down 
range sigma’s if necessary. 

The definite range law, as explained in the introduction, 
generally gives useful results and we feel it is quite 
appropriate for this problem. It may further be noted that 
this problem would not be solved were it not for this 
assump tion . 
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APPENDIX A 



A brief description of the integration technique employed. 

Again consider a 20 P*S*I. curve in 3-space. We will first 
enclose the surface in a "box", then choose, at random, a 
coordinate inside this box. If the coordinate is inside the 
surface of revolution we record the "height" of the function 
at that point; otherwise, we record a zero. The sum of these 
observations divided by the number of observations will give 
the "average height" of the function over the surface. The 
average height multiplied by the volume of the box will 
approximately be equal to the integral of the function con- 
cerned. As the number of trials gets large the approximation 
will, of course, be better. 

If we consider each coordinate chosen as a Bernoulli 



p = P r (coordinate is under surface) 

q = 1-p 

n = # of trials 

The standard deviation is (p 

In our case we chose n = 10,000 which yields an upper bound 
(p = q = h) on the deviation asl/% : \ : 10 ,000 “ *005 or a 
3 (f~ ’accuracy* of .015 



trial, the "worst" the variance could 




where 
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We chose : 



Yi, Y 2 , Y 3 = Uniform 


(0,1) random 


variables 


Xi = 2 c(Y x - 0 . 5 ) 


Xi : 


Uniform (- c, c ) 


X 2 = 2 c(Y 2 - 0 . 5 ) 


X 2 : 


Uniform (- c, c ) 


X 3 = b Y 3 


cn 

X 


Uniform ( 0 , b ) 


( b = 1000, c = 2000 


in our case) 





If the triple of random numbers (X^ , X2, X3) drawn satisfy 

2 2 2 

the inequality X-^ + X 9 < g (X^)the function is evaluated at 
the point (X]^ , X2 , X3) , otherwise a zero is recorded. 

Let 



2 2 2 

h(X-^) = integrand evaluated at (Xi , X2, X3) if X-^ + X2 — g (X3) 

. ^ otherwise 

we know that 



lim 

<Y\ OO 



1 £ 1, (x t ) E[A(x)J- f 



/ = • 



where ^(X) = 1 in our case. The required integral is therefore 

approximated by 2 - V* /v N 

be n h (X.) 

1 -a 



USE OF PROGRAM 

The following procedure is recommended for f inding^^^and PKILL: 

1 . Function X GOF must be rewritten to fit destruction curves 
employed . 

2 . Values of (p and(f) must be changed in the main program to 
suit the weapon employed. 
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3. Since only curves for a 1KT weapon are generally available, 
the problem must first be solved for this case using as a first 



approximation the value of j ^ ) found by the method of 
section IV and the parameters^ j. and^ must be multiplied by 

-V ) 

((f is not a function of the yield since it is inherent 
in the fusing mechanism). 



Conditional probability of destruction given the "best” 
trajectory when 



THETA 

A 

ALPHA 

^max 

XCOORD 

f 



= angle of re-entry 

= lower intercept of trajectory with destruction envelope 
= upper intercept of trajectory with destruction envelope 
= vertical coordinate of IBP 
* down range horizontal coordinate of IBP 
= down range distance from target the trajectory 
intercepts the ground 
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APPENDIX B 



PROGRAM PROB 
TYPE REAL MU 
DIMENSION MU{ 50) ,SI( 50) 

PRINT n 
1 1 FORMAT ( 1 HI ) 

PRINT 14 

14 FORM AT l 4X , 1 H 1 , 5X » 5HMU ,5X, 5HSI ,10HEXP0FXl 
1 » 1 OHP K ILL ) 

SI(l) = 670. $ MU ( 1 ) = 425. $1=1 
A=.(2.«3. 1415926535)*«1.5 
B=1 ./A 

152 THETA = 40. /57. 295779513 1 

153 SIGMA =200. $ F =SIGMA*SISMA $ G=l./F 

CC=TANF{ THETA) $ C = l./CC 

SIGMA1 = 500. 

D = S IGM A1 «S IGMA 1 
E = 1 ./D 

DO 777 LL = 1,70 
N = 0 

VS1 = VS2 = VS3 = 0. 

160 CONTINUE 

Z1 = RANF(-l) $ Z2 = R AN F (-2 ) $ Z3 = RANF (-3) 
XI = 2000. * ( Z 1 — . 5) $ X2 = 2300.MZ2-.5) $ X3 = 
Y 1 = X 1 *X1 $ Y2 = X2*X2 $ Y3 = X3*X3 
1023 W= XGOF ( X3 ) 



WN=W*W 

IFIY1+Y2.LE. WW ) 101,102 
102 V = VI = V2 = 0. $ GO TO 575 

10 1 CONTINUE 
H =S IGMA*D 
P=1 . /H$Q=B*P 
161 RR = Xl-SK I ) + X3*C 

R= RR*RR$S =X 2*X 2$ TT=X 3-MU { I )$T = TT*TT 
167 U=EXPF(-.5M ( E*( S+R) )+T*G) ) 

V = G*U $ VI = Xl*V $ V2 = X3*V 
575 CONTINUE 

VS1 = VS1+V $ VS2 = VS2+V1 $ VS3 = VS3+V2 
N = N+l 

IFtN.LT. 10000) 160,500 
500 XN=N 

P K ILL = 40000000C0.«VS1/XN 
El = 4000000000. «VS2/XN 
E3 = 4000000000. *VS3/XN 
PRINT 15, I, MU (I) »SI( I),E1,E3,PKILL 
15 FORMAT { 1X,I4,F10.5,F10.5,F10.5,F10.5,F10.5) 
I = 1 + 1 

HUH) = E3/PKILL 

SI(I) = (El+MUd )*PKILL*C)/ : >KILL 
777 CONTINUE 
END 



, 1 0 HEXP0FX3 



1 000. *Z3 



FUNCTION XGOF { Y ) 

1022 B 1 = 71.1 229667E Cl 
B2=-31 .7485590E-C1 
B3= 66.931 1997E-C3 
B4=-46. 251801 8E-C5 
B5= 1 4. 895 1 096 E-C7 
86= -2 4. 1298172 E- 10 
B7= 19.05551 0 1 E- 13 
B8 =-58. 4969877 E- 17 

XGOF = ({({((( B8*Y)+B7)*Y+B6)*Y+B5)*Y+B4)*Y+B3)*Y+B2)*Y+B1 

END 



APPENDIX C 



f. 



THETA 




A 


ALPHA 


Pmax 


XCOORD 


. 5.;- • 




J ' 

- l \ ! 

9 

10 

i 1 1 ‘ 1 


499 

424 

389 

369 

35l 

329 


.999 

.907 

.185 

.360 

.793 

.999 


•• '659.945 
660.099 
653.707 
646.386 
640.607 

636.991 


587’ 

531 

505 

494 

487 

479 


30.023 
60. 080 
46.210 
( 1.486) 
47.129 
77.037 


-6455° 

-4899 

-3865 

-3178 

-2712 

-2382 




1 1 
12 

13 

14 

J 5 

1 0 


952.323 

269.999 
235.552 
201.599 

169.999 


~ 6 _ 3y.204“ 
63^.545 
634.304 
633.962 
633.261 


— 45B 
454 
438 
420 
401 


90o03 

92.987 

91.381 

91.176 

95.059 


-2x30. 

-1923 

-1738 

-1564 

-1399 




17 

18 

19 

20 
•J \ 


iVl # J 

114.999 

89.067 

61.310 

29.999 


O 3 • c JL O 

631.038 

630.126 

629.931 

630.912 


36.9 

356 

346 

339 


16*25720"" 

111.540 

117.810 

118.230 

111.320 


“-T2*2 “ 

-1096 
— 964 
-849 
-753 




22 

23 

24 
_25 


w-j 

-00 

.•00 

-00 

QQ 


*vuu 

<.000 

*000 

.000 

..QQ0 


fi3T473— 
* 637.916 
644.412 
652.994 . 
663.553 


33* 

333 

334- 

336 

341 


567440 

84.139 

75.475 
80. 023 
102.240 


-070 
-6 20 
-583 
-562 
-560 ' 




27 

28 

29 

30 




075T8 59 
689.581 
704.321 
719.641 
* 735.102 


3*0 

352 

359 

366 

374 


137.9OO" 
1(56.230 
139.500 - 
-08. 566 
-78.982 


-562 

-576 

-596 

-620 

-644 




32 

33 

34 

35 




75”0 . 2o9 
764.841 
778.468 
790.967 
802.227 


381 

388 

395 

402 

408 


—>'6 . 0*3 
-103.780 
-117.540 
-131.330 
-144.830 


-670 

-687 

-704 

--716 

-725 




37 
; 38 

39 

40 




812. 226 :• 
.. 821.032 
828.780 
835.665 

841.917 


• 415 
421 
427 
432 
_ 437 


-I37.0YO 
-170.220 
-182.160 
-193.820 
-205.430. , 


"7779 

-730 . 

-728 
-724 
-719 




^1 

42 

43 

44 

41 




8477776 
853.488 
859.259 
865.257 
871. 596 


*42 
446 
450 
• 454. 
457 


— -21/72 20 
r 229. 370 ' 
-241.960 
-254.920 
-263.120 


4-7 1* 
-709 
-705 
-702 
-700 


. 


47 

48 
49. 
50 


\ 


/ 


8/6. 3 18 
885.401 
892.748 
900.205 
907.564 


“”460 

462 

463 

465 

466 


-281. 25O 
-293.940 
-305.720 
-316.090 
-324.550 


-699 

-700 

-700 

-700 

-700 
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APPENDIX C (Con't) 



THETA PROBABILITY OF DESTRUCTION 

SIGMA a 10 50 100 ’ 200 500 1000 

5 1.000 .890 .576 .311 .127 .064 



6 

7 

8 

9 

10 


.981 .760 .443 .186 .094 

.992 .814 ' .492 .209 .105 

.994 .834 1511 .218. .110 

.996 .851 i530 .227 .115 

- - - 1948 -87 4 :W7 • .741 


11 

12 

% 

15 


l.C 


f99 .904 .595 .261 .132 

)00 .932 .638 .285 .145 

.954 .681 .310 .158 

.969 .720 .335 * .171 

.979 -777 -777 -187 


lb 

17 

18 

19 

20 




.986 .780 .376 Tl94 
.990 .803 .394 .204 
.993 .824 .412 *.213 
.99 6 .845 .430 .224 

.997 .867 .477 .7 76 


21 

22 

u 

— . - 25 




.999 .889 .477 .250 

.999 .909 .502 .265 

1.000 .923 .521 .277 

1.000 .925 .524 .278 

.999 .903 .497 .7*0 


^0 

27 

28 

29 

30 




.9 

.9 

1.0 


99 .909 • .501 .265 

99 ‘ .915 *510 .270 

00 .922 .519 .275 

.928 .528 .281 

- - /934 ...... . 778 .787 ... 


31 

32 

S 

' 35 






.939 .5^7 .292 

.944 .556 • .298 

.948 .564 .303 

- .952 .571 .308 

_ .977. . 778 717. 


363 

37 

33 

2? 






.958 .583 .315 

.960 .588. .319 

.962 .593 .321 

.963 • .5977 .324 

_ . *90^ _ . .. *996 


41 • 

42 

tl 






.966 .603 .328 

.967 .607 .330 

.968 .610 .323 

.969 .613 .335 

"971 .. .7,1 7 .777 ... 


4o 

4? . 

48 

49 

50 V 


1 

' \ 


' ^ 


• .972 .620 .339 

.973 .624 .342 

.974 .628 .345 

.976 .632 .347 

' .977 .636 .350 
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APPENDIX C (don't) 



THETA 


A 


ALPHA 


/^max 


51 


00.000 


909.663 


466 


52 




911.767 


. 467 


53 




921.521 


. 467 


54 




930.684 


467 


55 




932.232 


468 



3T 

57 

58 

P 

6 0 



XCOORD 

-330 

-333 

-33$ 

-331 

-324 



f 

6 95 
689 
681 
669 
653 



ST" 

62 

63 

64 

65 



933.000 


468 


-315 


633 


933.000 


468 


-303 


611 




467 


-289 


585 




467 


-273 


• 556 




467- 


-256 


526 



467 


-239 


496 


466 


-223 


466 


465 


-207 


438 


465 


-194 


412 


464 


-I83 


389 



66 


464 


-174 


371 


67 


464 


-169 


■ 358 


68 


464 


-165 


345 


69 


465 


-164 . 


338 


70 


467 


-164 • 


333 



71 

72 

73 

74 

75 



75" 

77 

78 

79 

80 



470 


-164 


330 


474 

479 


z\U 


326 

322 


484 


-160 


314 


490 


-154 


301 



494- 


-144 


283 


495 


-130 


259 


493 

484 


-113 

-937 


230 

196 


468 


-732 


160 



00 00 00 00 00 
Vn N> m 


“ 446 _ 
401 
350 
288 
221 


- 269 

— 24l 
-297 


126 

986' 

801 

756 

865 


S5 


' 158 


-430 


111 


87 


116 


-588 


142 


88 


124 


-682 


161 


89 


221 


-538 


136 


90 


468 


-105 


165 
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APPENDIX C (Con't) 



TKETA 



SIGMA 



PROBABILITY OF DESTRUCTION 
10 50 100 200 500 



1000 



52 j 

It 

55 






. . 9 ^ 

.9' 
• 91 


79 .643 .355 

79 .646 .357 / 

30 .648 .358 • 

.649 *359 


55 

57 

58 

I 9 
6 0 


• 




- 


. 


4-9 

i 


. 


61 

62 
o3 
64 

o 5 






1 






• 


oo 

67 

68 ■ 

69 

70 


• 




1 








71 

72 

73 

74 

75 








1 




- 


76TT" 

77 

78 

79 

80 














6l 

82 

83 

84 

85 










* / t 


; * 


85 

87 ' 

88 . 

89 

90 . \ 


' \ 


r. \ 


i 

i 

i 

' 


- ■ i 


( ' \ 
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